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Abstract — We consider a decentralized bidirectional control of 
a platoon of A'^ identical vehicles moving in a straight line. The 
control objective is for each vehicle to maintain a constant velocity 
and inter-vehicular separation using only the local information 
from itself and its two nearest neighbors. Each vehicle is modeled 
as a double integrator. To aid the analysis, we use continuous 
approximation to derive a partial differential equation (PDE) 



approximation of the discrete platoon dynamics. The PDE model 



is used to explain the progressive loss of closed-loop stability with 
O increasing number of vehicles, and to devise ways to combat this 
loss of stability. 

If every vehicle uses the same controller, we show that the least 
stable closed-loop eigenvalue approaches zero as O(-j^) in the limit 
of a large number (A'^) of vehicles. We then show how to ameliorate 
I Xhis loss of stability by small amounts of "mistuning", i.e., changing 
the controller gains from their nominal values. We prove that with 
arbitrary small amounts of mistuning, the asymptotic behavior of 
the least stable closed loop eigenvalue can be improved to O( j^). 
All the conclusions drawn from analysis of the PDE model are 
corroborated via numerical calculations of the state-space platoon 
2 model. 
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\ I. Introduction 

We consider the problem of controlling a one-dimensional 
platoon of N identical vehicles where the individual vehicles 
move at a constant pre-specified velocity Vd with an inter- 
vehicular spacing of A. Figure [Tla) illustrates the situation 
schematically. This problem is relevant to automated highway 
systems (AHS) because a controlled vehicular platoon with a 
constant but small inter-vehicular distance can help improve 
the capacity (measured in vehicles/lane/hour, as in [JJ) of a 
highway JU . Due to this, the platoon control problem has been 
extensively studied |31, [4|, |51, [l], [6|, |7|. The dynamic and 
control issues in the platoon problem are also relevant to a 
general class of formation control problems including aerial 
V^ehicles, satelUtes etc. El, ID. 

Several approaches to the platoon control problem have been 
considered in the literature. These approaches fall into two broad 
categories depending on the information architecture available 
to the control algorithm(s): centralized and decentralized. We 
call an architecture decentralized if the control action at any 
individual vehicle is computed based upon measurements ob- 
tained by on-board sensors and possibly wireless communication 
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with a limited number of its neighbors. We call all other ar- 
chitectures centralized. Decentralized architectures investigated 
in the literature include the predecessor-following |[T|, ifTOl . 
1111 and the bidirectional schemes |12l, f?!, |T3l, |T4l, fBI. 
In the predecessor-following architecture, the control action at 
an individual vehicle depends only on the spacing error with 
the predecessor, i.e., the vehicle immediately ahead of it. In 
the bidirectional architecture, the control action depends upon 
relative position measurements from both the predecessor and 
the follower On the other hand, in a centralized architecture 
measurements from all the vehicles may have to be continually 
transmitted to a central controller or to all the vehicles. The 
optimal QR designs of f4l, f6l typically lead to centralized 
architectures. Predecessor and Leader follower control schemes 
(see lfT6l . ifTTl and references therein), which require global 
information from the first vehicle in the platoon are also ex- 
amples of the centralized architecture. The high communication 
overhead in a centralized architecture makes it less attractive for 
platoons with a large number of vehicles. Additionally, with any 
centralized scheme, the closed loop system becomes sensitive 
to communication delays that are unavoidable with wireless 
communication 1 1 8 1 . 

The focus of this paper is on a decentralized bidirectional 
control architecture: the control action at an individual vehicle 
depends upon its own velocity and the relative position errors 
between itself and its predecessor and its follower vehicles. The 
decentralized bidirectional control architecture is advantageous 
because, apart from its simplicity and modularity, it does not 
require continual inter-vehicular communication. Measurements 
needed for the control can be obtained by on-board sensors 
alone. Each vehicle is modeled as a double integrator. A double 
integrator model is common in the platoon control literature 
since the velocity dependent drag and other non-linear terms 
can usually be eliminated by feedback linearization ||T|, ifTOl . 
The control objective is to maintain a constant inter-vehicular 
spacing. 

In spite of the advantages over centralized control, there 
are a number of challenges in the decentralized control of a 
platoon, especially when the number of vehicles, N, is large. 
First, the least stable closed-loop eigenvalue approaches zero 
as the number of vehicles increases |fT9l . Among decentralized 
schemes, one particularly important special case is the so- 
called symmetric bidirectional control, where all vehicles use 
identical controllers that are furthermore symmetric with respect 
to the predecessor and the follower position errors. In this 
case, the least stable closed loop eigenvalue approaches as 
0{^) with a symmetric bidirectional control and this behavior 
is independent of the choice of controller gains |19|. This 
progressive loss of closed-loop damping causes the closed loop 
performance of the platoon to become arbitrarily sluggish as 
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the number of vehicles increases. It is interesting to note that 
the 0( decay of the least stable eigenvalue occurs with the 
centralized LQR control as well 16\. 

The second challenge with decentralized control is that the 
sensitivity of the closed loop to external disturbances increases 
with increasing N. With predecessor following control, distur- 
bances acting on the vehicles cause large inter-vehicular spacing 
errors Ol, |[T1, ll20l The seminal work of |20| on string instability 
was partly inspired by this issue. It was shown in Q that 
sensitivity to disturbances with predecessor following control is 
independent of the choice of the controller Similar controller- 
independent sensitivity to disturbances is also exhibited by the 
symmetric bidirectional architecture Q, ifTSl . 11211 . In |22|, it 
was shown that symmetric architectures have similarly poor 
sensitivity even when every vehicle uses information from more 
than two neighbors, as long as the number of neighbors is no 
more than 0{N'^^^). 

Third, there is a lack of design methods for decentralized 
architectures. For N vehicles, in general, N distinct controllers 
need to be designed, for which few control design methods 
exist. This has led to the examination of only the symmetric 
control among bidirectional architectures Q, |fT3l , ll22l . Some 
symmetry aided simplifications are possible for analysis and 
design in this case. 

In summary, while issues such as stability and sensitivity to 
disturbances become critical as the platoon size increases, a lack 
of analysis and control design tools in decentraUzed settings 
makes it difficult to address these issues. 

In this paper we present a novel analysis and design method 
for a decentralized bidirectional control architecture that amelio- 
rates the progressive loss of closed loop stability with increasing 
number of vehicles. There are three contributions of this work 
that are summarized below. 

First, we derive a partial differential equation (PDE) based 
continuous approximation of the (spatially) discrete platoon 
dynamics. Just as a PDE can be discretized using a finite 
difference approximation, we carry out a reverse procedure: 
spatial difference terms in the discrete model are approximated 
by spatial derivatives. The resulting PDE yields the original set 
of ordinary differential equations upon discretization. 

Two, we use the PDE model to derive a controller independent 
conclusion on stability with symmetric bi-directional architec- 
ture. In particular, the behavior of the least stable eigenvalue 
of the discrete platoon dynamics is predicted by analyzing the 
eigenvalues of the PDE. We show that the least stable closed- 
loop eigenvalue approaches zero as 0{-^). This prediction is 
confirmed by numerical evaluation of eigenvalues for both the 
PDE and the discrete platoon model. The real part of the least 
stable eigenvalue of the closed loop is taken as a measure of 
stability margin. 

The third and the main contribution of the paper is a 
mistuning-based control design that leads to significant improve- 
ment in the closed loop stability margin over the symmetric 
case. The biggest advantage of using a PDF-based analysis is 
that the PDE reveals, better than the state-space model does. 



the mechanism of loss of stability and suggests a mistuning- 
based approach to ameliorate it. In particular, analysis of the 
PDE shows that forward-backward asymmetry in the control 
gains is beneficial. The asymmetry refers to the assignment of 
controller gains such that a vehicle utilizes information from the 
preceding and following vehicles differently. Our main results. 
Corollary |2] and Corollary |3] give control gains that achieve 
the best improvement in closed-loop stability by exploiting this 
asymmetry. In particular, we show that an arbitrarily small per- 
turbation (asymmetry) in the controller gains from their values 
in the symmetric bidirectional case can result in the least stable 
eigenvalue approaching only as (as opposed to 0{^) 

in the symmetric bidirectional case). Numerical computations 
of eigenvalues of the state-space model of the platoon is used 
to confirm these predictions. Mistuning based approaches have 
been used for stability augmentation in many applications; 
see 11231 . Il24l . ||25| . Il26l for some recent references. Our 
paper is the first to consider such approaches in the context 
of decentralized control design. 

Although the PDE model is derived under the assumption of 
large N, in practice the predictions of the PDE model match 
those of the state-space model accurately even for small values 
of N . Similarly, the benefits of mistuning are significant even 
for small values of N (see Section IVII i. 

In addition to the stability margin improvements, the mis- 
tuning design reduces the closed loop's sensitivity to external 
disturbances as well. In bidirectional architectures, the 
norm of the transfer function from the external disturbances 
to the spacing errors is used as a measure of sensitivity to 
disturbances [7J. Numerical computation of the Hoc norm of 
this transfer function shows that mistuning design also reduces 
sensitivity to disturbances significantly (see Section FVI-Db . 

We briefly note that there is an extensive literature on mod- 
eling traffic dynamics using PDEs; see the seminal paper of 
Lighthill and Whitham |27| for an early reference, the paper 
of Helbing |28| and references therein for a survey of major 
approaches, and the papers of [29| and |30| for control-oriented 
modeling. In spite of apparent similarities, our approach is 
quite different from the existing approaches. PDE models of 
traffic dynamics typically start with continuity and momentum 
equations (28^. Moreover, one requires a model of human 
behavior to determine an appropriate form of the external force 
in the momentum equation. This difficulty frequently leads to 
the introduction of terms in the PDE that are determined by 
fitting data; see [28, Section III-D] for a thorough discussion of 
such approximations used in various continuum traffic models. 
In contrast, we approximate the closed loop dynamic equations 
by continuous functions of space (and time) that are inspired by 
finite-difference discretization of PDEs. Ad-hoc approximations 
of human behavior is not needed. Moreover, the original dy- 
namics can be recovered by discretizing the derived PDE, which 
provides further evidence of consistency between the (spatially) 
discrete and continuous models. 

We also note that macroscopic models of traffic flow models 
have been used for designing control laws for a complete 
automated highway system (AHS) with lane changing, merging, 
etc. in addition to a platoon in one lane (see OOl . lISTI and 
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(a) A platoon with fictitious lead and follow vehicles. 




(b) Same platoon in y coordinates. 
Fig. 1. A platoon with vehicles moving in one dimension. 



references therein). The PDE model derived in the paper is not 
applicable to a complete AHS, but only to a single platoon. 

The rest of the paper is organized as follows: section HH states 
the platoon problem in formal terms by describing a state-space 
model of the closed loop platoon dynamics; section |III] then 
describes the derivation of the PDE model from the state space 
model. In section HV] the PDE is analyzed to explain the loss of 
stability with N, and section|V]describes how to ameliorate such 
loss of stability by mistuning. Section IV-CI reports simulation 
results that show the benefit of mistuning in time-domain. In 
Section |Vll we comment on various aspects of the proposed 
mistuning design. 

II. Closed loop dynamics with bidirectional 

CONTROL 

Consider a platoon of N identical vehicles moving in a 
straight line as shown schematically in Figure[Tla). Let Zi{t) and 
Vi{t) :— Zi{t) denote the position and the velocity, respectively, 
of the i* vehicle for i = 1, 2, . . . , A^. Each vehicle is modeled 
as a double integrator: 



(1) 



where Ui is the control (engine torque) applied on the i* vehicle. 
Such a model arises after the velocity dependent drag and 
other non-linear terms have been eliminated by using feedback 
linearization lH], ifTol . 

The control objective is to maintain a constant inter-vehicular 
distance A and a constant velocity Vd for every vehicle. Every 
vehicle is assumed to know the desired spacing A and the 
desired velocity Vd- The control architecture is required to 
be decentralized, so that every vehicle uses locally available 
measurements. We assume that the error between the position 
(as well as velocity) of a vehicle and its desired value is small, 
so that analysis of the platoon dynamics with a linear vehicle 
model and a linear control law is justified. 

In this paper, we assume a bi-directional control architecture 
for individual vehicles in the platoon (except the first and the 
last vehicles). For the first and the last vehicles, we consider two 
types of control architectures (termed as scenarios I and II) as 
tabulated in Table U In scenario I, we introduce (after ([6], Q) 



Scenario 


Length L 


Leader 


Follower 


I 


(iV + l)A 


Co = 


Cjv+i = 


II 


7VA 


■Do = 





TABLE I 
The two scenarios. 



a fictitious lead vehicle and a fictitious follow vehicle, indexed 
as and + 1 respectively. Their behavior is specified by 
imposing a constant velocity trajectories as Zo(i) = Vdt and 
Zjv+i — Vdt — {N + 1) A. In scenario II, only a fictitious lead 
vehicle with index i — with Zq (t) = Vdt is introduced. For the 
last vehicle in the platoon in scenario II, there is no follower 
vehicle and it uses information only from its predecessor to 
maintain a constant gap. 

Consistent with the decentralized bidirectional linear control 
architecture, the control Ui for the i* vehicle is assumed to 
depend only on 1) its velocity error Vi — Vd, and 2) the relative 
position errors between itself and its immediate neighbors. That 
is. 



kf\z. 



Z.^+1 - A) 
~h{V~Vd). 



(2) 



where fc| ^ bt are positive constants. The first two terms are used 
to compensate for any deviation away from nominal position 
with the predecessor (front) and the follower (back) vehicles 
respectively. The superscripts and '^''^ correspond to front 
and back, respectively. The third term is used to obtain a zero 
steady-state error in velocity. In principle, relative velocity errors 
between neighboring vehicles can also be incorporated into the 
control, but we do not examine this situation here. Since Vd 
and A are known to every vehicle, the relative errors used in 
the control law, including the velocity error, can be obtained in 
practice by on-board devices such as radars, GPS, and speed 
sensors. 

The control law (O represents state feedback with local (near- 
est neighbor) information. Analysis of this controller structure 
is relevant even if there are additional dynamic elements in the 
controller There are several reasons for this. First, a dynamic 
controller cannot have a zero at the origin. It will result in a pole- 
zero cancellation causing the steady-state errors to grow without 
bound as A^ increases lfT3l . Second, a dynamic controller cannot 
have an integrator either For if it does, the closed-loop platoon 
dynamics become unstable for a sufficiently large values of 
A^ |13|. As a result, any allowable dynamic compensator must 
essentially act as a static gain at low frequencies. The results 
of lfT3l indicate that the principal challenge in controlling large 
platoons arises due to the presence of a double integrator with its 
unbounded gain at low frequencies. Hence, the limitation and 
its amelioration discussed here with the local state feedback 
structure of (O is also relevant to the case where additional 
dynamic elements appear in the control. 

To facilitate analysis, we consider a coordinate change 



Vi = 27r( 



Z,{t)-Vdt + L^ 



2n 



L ' 



(3) 
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where L denotes the desired platoon length, which equals {N + 
1)A in scenario I and A^A in scenario II. Figure [TJb) depicts 
the schematic of the platoon in the new coordinates. The scaling 
ensures that yQ{t) = 2tt, yi{t) G [0, 27r], and yM+i{t) = 
iyN{t) — 0) in scenario I (II). Here, we have implicitly assumed 
that deviations of the vehicle positions and velocities from their 
desired values are small. 

In the scaled coordinate, the dynamics of the vehicle are 
described by 

Vi = Ui, (4) 
where Ui := 2iTUi/L. The desired spacing and velocities are 
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0, 



(5) 



L/27r' L/2t: 
and the desired position of the vehicle is 

yf{t) = 27T- iS. (6) 
The position and velocity errors for the i* vehicle are given by: 

yi(t) = Uiit) - yfit), Vi = Vi -Vd = Vi, and 

yi = Vi- 



(7) 



We note that vq — vn+i = for the fictitious lead and 
follow vehicles. In the scaled coordinates, the decentralized 
bidirectional control law ^ is equivalent to the following 



,(/) 



Vi+i - 5) 



{Vi-i - Vi) - kf {jji - yi+i) - h,v. 



(8) 



It follows from (01) and ^ that the closed loop dynamics of the 
vehicle in the y-coordinate is 



Vi + b,y, = k\^\y^_i - iji) 



kf\y',-yt+i)- 



(9) 



To describe the closed-loop dynamics of the whole platoon, we 
define 



,yN\ 



,vn\ 



For scenario I with fictitious lead and follow vehicles, the control 
law (O yields the following closed loop dynamics. 



(/) 







B 



(10) 



where 

diag(fcf\A;f 



= diag{k[^\k!/\ 

Jb) 



K 



(b) 



, ky), B = diag(5i,62, . . . ,6jv), and 



M = 



1-10 
1-1 





1 -1 

1 



For scenario II with a fictitious lead vehicle and no follow 
vehicle, the closed loop dynamics are 







I ■ 

-B 



(11) 



where k[P = KY', K\1' ^ diag(fci°\fe 



(/) 1^(6) 



Sb) Ab) 



, A:^Li,0), and 



1 -1 
1-1 





1 -1 

. 



Our goal is to understand the behavior of the closed loop 
stability margin with increasing N and to devise ways to 
improve it by appropriately choosing the controller gains. While 
in principle this can be done by analyzing the eigenvalues of 
the matrix A^^p (scenario I) and of Al (scenario II), we 
take an alternate route. For large values of N, we approximate 
the dynamics of the discrete platoon by a partial differential 
equation (PDF) which is used for analysis and control design. 

III. PDF MODEL OF PLATOON CLOSED LOOP DYNAMICS 

In this section, we develop a continuous PDF approximation 
of the (spatially) discrete platoon dynamics. The PDF is derived 
with respect to a scaled spatial coordinate x E [0, 2tt]. We recall 
that in Section |II] the scaled location of the vehicle (denoted 
as yi) was defined with respect to such a coordinate system. 
In effect, the two symbols x and y correspond to the same 
coordinate representation but are used here to distinguish the 
continuous and discrete formulations. As in the discrete case, 
the platoon always occupies a length of 2tt irrespective of N. 

A. PDE derivation 

The starting point is a continuous approximation: 

v{x,t) := Vi{t) at X = yi 
^ t) = i^(t). (from ©) 

Similai-ly, fcf^) (a;), fc('') {x) are used to denote continuous 
approximations of discrete gains hi,k\^\kf'^ respectively. We 
will construct a PDF approximation of discrete dynamics in 
terms of these continuous approximations. To do so, it is 
convenient to first differentiate (|9]l with respect to time, 

Vi + hiVi = k\^\vi-i - Vi) - kf\vi - Vi+i). (12) 

We recast this equation 



— where 



It follows that 

Vi + biVi = ^kl'\vi-i - Vi+i) + ^k'^^\vi-i - 2vi + Vi+i) 
1 (_) Vi^i - Ui+i 1 (+) Vi-i - 2{i, + Vi+i 



— k 



26 2pl 



,2 '^i 



51 



where 



PO 



1 _ N 
'6~ 2^' 



(14) 
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Pq has the physical interpretation of the mean density (vehicles 
per unit length). Now, we make a finite-difference approximation 
of derivatives 

d 



Vi-l 



26 

2vi + Vi+i 



dx 



v(x, t) 



x=yi 



where we recall that v{x,t) is a continuous approximation of 
the vehicle velocities {vi{t) — v{yi,t) etc). Denoting k'^-^Ux) 
and k'^^\x) as continuous approximations of fcl^-* and ' 
respectively, the discrete model is written as: 



9*2 



v{x,t) 



d 

Kx)-Q-^v{x,t) 



1 

Po L 



k^-\x)—v{x,t) 



1 

Wo 



g2 

k^+\x)^v{x,t) 



Hence, we arrive at the partial differential equation (PDE) as a 
model of the discrete platoon dynamics: 



d 

+ b{x)—]vix,t) = 



d 



-k(-Hx)^ 

Po ox 



v{x,t) (15) 



In the remainder of this paper, we assume that fc(+)(.T) > 0. 
Using dTsT l. the continuous counterparts of the front and the 
back gains are given by 



k^f\x) = i (/c(+)(a;) + fc(-)(x)) , 



(16) 



k^fHy^ 



so that the gain values fc,- can be obtained as k'^'^'' 
and /cf = k'^^'^y^). It can be readily verified that one recovers 
the system of ordinary differential equations (( fT2] i for i — 
1 , . . . , A^) by discretizing the PDE ( fTSl ) using a finite difference 
scheme on the interval [0, 2tt] with a discretization 6 between 
discrete points. 

The boundary conditions for the PDE ( fTsT i depend upon the 
dynamics of the first and the last vehicles in the platoon. For 
scenario I with a constant velocity fictitious lead and follow 
vehicles, the appropriate boundary conditions are of the Dirichlet 
type on both ends: 



v{0,t) = v{27T,t)^0, V<e[0,oo). 



(17) 



For scenario II with the only a fictitious lead vehicle, the 
appropriate boundary conditions are of Neumann-Dirichlet type: 

dv 



dx 



(0,i) = w(27r,t) = 0. Vte[0,oo) 



(18) 



We refer the reader to Appendix II-AI for a discussion on well- 
posedness of the solutions to ( fTSl l. It is shown in Appendix II-AI 
that a solution exists in a weak sense when fc(+\ k^^\ ^lET ^ 
L-([0,27r]). 

Equation ( fTSl ) describes spatio-temporal evolution of small 
velocity perturbations in a platoon. It is worthwhile to note that 
the PDE model is a hyperbolic equation. Without the two first 



platoon 
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pde 
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-0,5 
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(a) Scenario I ( Dirichlet- 
Dirichlet) 



(b) Scenario II ( Neumann- 
Dirichlet) 



Fig. 2. Comparison of closed loop eigenvalues of the platoon dynamics and 
the eigenvalues of the corresponding PDE 1191 for the two different scenarios: 
(a) platoon with fictitious lead and follow vehicles, and con'espondingly the 
PDE )19t with Dirichlet boundary conditions, (b) platoon with fictitious lead 
vehicle, and correspondingly the PDE il9\ with Neumann-Dirichlet boundary 
conditions. For ease of comparison, only a few of the eigenvalues are shown. 

f f) 

Both plots are for N = 25 vehicles; the controller parameters are fc" ' = 
k^''^ = 1 and bi = 0.5 for i = 1,2, . . . , Af, and for the PDE k'-f^x) = 
k(*'){x) = 1 and b(x) = 0.5. 



order terms (i.e., for b{x) = k'--\x) = 0), the PDE is a standard 
wave equation with spatially inhomogeneous values of wave 
speed. The term ^fc*- ■'(^)ff is an advection term, and b{x)^ 
is a damping term. The hyperbolic nature of the PDE model 
means that a perturbation originating, say, in the middle of a 
long platoon will propagate both upstream and downstream with 
finite speed. The two first order terms serve to modify aspects 
of this propagation. The damping term causes a perturbation to 
damp out in time. The advection term serves to create possible 
asymmetries in upstream versus downstream propagation. 



B. Eigenvalue comparison 

For preliminary comparison of the PDE obtained above with 
the state-space model of the closed loop platoon dynamics, we 
consider the simplest case where the position control gains are 
constant for every vehicle, i.e., k'^^\x) = ^(''^(x) = fco and 
b{x) = bo. In such a case k^^\x) = 0, fc(+'(a;) = 2ko and the 
PDE (fTSl l simplifies to 



bo 



d_ _ ko \ 
dt Pq dx^ J 



0, 



(19) 



which is a damped wave equation with a wave speed of 
The wave equation is consistent with the physical intuition that a 
symmetric bidirectional control architecture causes a disturbance 
to propagate equally in both directions. 

Figure |2] compares the closed loop eigenvalues of a discrete 
platoon with iV = 25 vehicles and the PDE (fT9] l. The eigenval- 
ues of the platoon are obtained by numerically evaluating the 
eigenvalues of the matrices A^-p and Al (defined in ( fTOb and 
(fTTTl). The eigenvalues of the PDE are computed numerically 
after using a Galerkin method with Fourier basis |32|. The figure 
shows that the two sets of eigenvalues are in excellent match. In 
particular, the least stable eigenvalues are well-captured by the 
PDE. Additional comparison appears in the following sections, 
where we present the results for analysis and control design. 
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boundai'y condition 


eigenvalue A; 


eigenfunction 


I 


V{0) = »?(27r) = 
(Dirichlet - Dirichlet) 


(2 
4 


sin(^) 


1 = 1,2,... 


|H(0) = r,(2^)=0 
(Neumann - Diiichlet) 


(2l-lf 
16 


, (2l-l)x X 
cos( 4 ) 


1 = 1,2,... 



boundary condition 



Diriclilet-Dirichlet 



Neumann-Dirichlet 



for I << Ic Ic 



4i)o iV^ 



+ 0(7^) 



TABLE III 



TABLE II 

The eigen-solutions for the Laplacian operator with two 
different boundary conditions. 



IV. Analysis of the symmetric bidirectional case 

This section is concerned with asymptotic formulas for 
stabiUty margin (least stable eigenvalue) for the symmetric 
bidirectional architecture with symmetric and constant control 
gains: k'^^^x) = k'-''\x) = fco and b{x) = bo- The analysis is 
carried out with the aid of the associated PDE model: 



where x e [0, 27r] and 



fco 

Pi 



0, 



(20) 



(21) 



(22) 



is the wave speed. The closed-loop eigenvalues of the PDE 
require consideration of the eigenvalue problem 

dx"^ 

where 77 is an eigenfunction that satisfies appropriate boundary 
conditions: ( fTTb for scenario I and ( fTSb for scenario II. The 
eigensolutions to the eigenvalue problem ( |23] ) for the two 
scenarios are given in Table [III The eigenfunctions in either 
scenario provide a basis of L^{[Q, 27r]). 

After taking a Laplace transform, the eigenvalues of the PDE 
model ( I20I 1 are obtained as roots of the characteristic equation 



60S 



0, 



(23) 



where A satisfies ( l22b . Using Table [III these roots are easily 
evaluated. For instance, the ^* eigenvalue of the PDE ( |20] i with 
Dirichlet boundary conditions is given by 



-bo ± V^o - «o^' 



(24) 



where 1 = 1,2,.... The real part of the eigenvalue depends upon 
the discriminant D{l,N) {bg — UqP), where the wave speed 
oo depends both on control gain fco and number of vehicles 
N (see (I2TI)). For a fixed control gain, there are two cases to 
consider: 

1) If D{l,N) < 0, the roots are complex with the real 



part given by — ^, 
2) If D{1, N) > 0, the roots sf are real with 



-bo- 



In the former case, the damping is determined by the velocity 
feedback term 60^, while in the latter case one eigenvalue (s^) 
gains damping at the expense of the other (sf) which looses 
damping. When sf^ are real, the eigenvalue is closer to the 



origin than s^; so we call the less-stable eigenvalue. 
The following lemma gives the asymptotic formula for this 
eigenvalue in the limit of large N. 

Lemma 1: Consider the eigenvalue problem for the linear 
PDE ( I20I 1 with boundary conditions ( [TtI i and (fTSl l. corresponding 
to scenarios I and II respectively. The /* less-stable eigenvalue 
approaches as 0{1/N^) in the limit as ^ cx). The 



asymptotic formulas appear in Table | 



□ 



Proof of Lemma |7] We first consider scenario I with Dirichlet 
boundary conditions (fTTI l. Using (l24l l and (ISTT i. 



2s± = -bo ± boll 



-bo ±bo{l 



27r2fco f 



o{ 



for alP/bl 
numbers 



<< 1. The asymptotic formula holds for wave 

boN 



ho_ 

a-o 



27^^/fci 



(25) 







and in particular for each / as A^ — ^ cxd. The proof for the 
scenario II with Neumann-Dmchlet boundary conditions (fTSl) 
follows similarly. 

The stability margin of the platoon can be measured by the 
real part of s'^ , the least stable eigenvalue. 

Corollary 1: Consider the eigenvalue problem for the linear 
PDE ( I20I ) with boundary conditions ( fTTI l and ( fTST l, corresponding 
to scenarios I and II respectively. The least stable eigenvalue, 
denoted by s'l, satisfies 



TT^fco 



1 



60 
TT^fco 



1 



Abo 



(Dirichlet-Dirichlet) 
(Neumann-Dirichlet) 



as A^ 



(26) 



(27) 



□ 



The result shows that the least stable eigenvalue of the closed 
loop platoon decays as with symmetric bidirectional control. 

We now present numerical computations that corroborates 
this PDF-based analysis. Figure |3] plots as a function of A^ the 
least stable eigenvalue of the PDE and of the state-space model 
of the platoon, as well as the prediction from the asymptotic 
formula. The eigenvalues for the discrete platoon are obtained 
by numerically evaluating the eigenvalues of the matrices ^l-f 
and Al (see (fTOb and ( fTTT )) with constant control gains k^'^'' = 



fco 



1 and hi 



bo 



0.5 for i = 1, 



,N. The 
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PDE (Holi, im 
platoon, L-F 
eq. l|26ll 



• PDE (EOll, N-D 

platoon, L 
-eq. l l27l l 




Fig. 3. Comparison of the least stable eigenvalue of the closed loop platoon 
dynamics and that predicted by Corollary[T]with symmetric bidirectional control. 
There are three plots each for scenarios 1 and II (corresponding legends are 
boxed together), and those three should be compared with one another. In 
the plot legends, "D-D" stands for "Dirichlet-Dirichlet", "N-D" for "Neumann- 
Dirichlet", "L-F" for fictitious leader-follower, and "L" for fictitious leader. The 
plot for "PDE )20) , D-D" should be compared with "platoon, L-F" since they 
both correspond to scenario L Similarly, "PDE i20\ . N-D" and "platoon, L" 
correspond to scenario II. Note that the predictions (26) and 1271 are valid for 
1 << Ic (defined in (25)), which in this case means for N >> 12. 



Re 



Im 



Re 



Im 



(a) Eigenvalues move toward 
zero with increasing N. 



(b) Mistuning "exchanges" 
stability between sf and s7 . 



Fig. 4. A schematic explaining the loss of stability as N increases and how 
mistuning ameliorates this loss. 



comparison shows that the PDE analysis accurately predicts the 
eigenvalue of the state-space model of the platoon dynamics. 

Figure Ufa) graphically illustrates the destabilization by de- 
picting the movement of eigenvalues sf as N increases. For 
sufficiently small values of N, the discriminant D{1,N) is 
negative and the eigenvalue sf are complex. The real part of 
the eigenvalue depends only on the value of Bq. At a critical 
value of = iVc := the discriminant becomes zero, 

bo 

sf = and the eigenvalues collide on the real axis. For values 
of N > Nc and in particular as iV ^ oo, the eigenvalue sf 
asymptotes to while staying real, and sj~ asymptotes to —b. 



Their cumulative damping, as reflected in the sum sf 
—bo, is conserved. In other words, Sj* 
expense of Sj^. 



is destabilized at the 



Remark 1: The preceding analysis shows that the loss of 
stability experienced with a symmetric bidirectional architecture 
is controller independent. The least stable eigenvalue approaches 
as 0{1/N'^) irrespective of the values of the gains fco and bo, 
as long as they are fixed constants independent of N. Corollary [T] 
also implies that for the least stable eigenvalue to be uniformly 
bounded away from 0, one has to increase the control gain ko as 
iV^. In [61, the same conclusion was reached for the least stable 
eigenvalue with LQR control of a platoon on a circle. LQR 
control typically leads to a centralized architecture, whereas 
symmetric bidirectional control is decentralized. It is interesting 
to note that the least stable eigenvalue behaves similarly in these 
distinct architectures. □ 



V. Reducing loss of stability by mistuning 

In this section, we examine the problem of designing the 
control gain functions k^^\x) , k^'^^x) so as to ameliorate the 
loss of stability margin with increasing N that was seen in the 
previous sections when k^^\x) = fc'''' = ko- Specifically, we 
consider the eigenvalue problem for the PDE ( fTSI ) where the 
control gains are changed slightly (mistuned) from their values 
in the symmetric bidirectional case in order to minimize the 
least-stable eigenvalue sf. With symmetric bidirectional control, 
one obtains an 0{^) estimate for the least stable eigenvalue 
because the coefficient of term in PDE ( fTSl l is 0{^) and 
the coefficient of S- term is 0. Any asymmetry between the 
forward and the backward gains will lead to non-zero k'-^' (x) 
and a presence of term as coefficient of By a judicious 

choice of asymmetry, there is thus a potential to improve the 
stability margin from 0{^) to 0{jj). 

We begin by considering the forward and backward position 
feedback gain profiles: 

k^f'ix) = fco + efc^^'^"''*^(a;), 
fc('')(x) = fco+efc^^'^"''*^(x), 

where e > is a small parameter signifying the amount of 
mistuning and k^^'P^'^''^^ (x), fc^'''^"''*^ (x) are functions defined 
over the interval [0, 2tt] that capture perturbation from the 
nominal value fco. Define 

k,{x) := fc(^'P"''*)(a;) + fc(^'P"''*)(a;), 
kra{x) := fc(^'f"'-*)(a;) - fc(''^f"'-*)(2;), 

so that from (fTSI l. 

fc(+)(a;) = 2fco + efcs(a;), fc(")(a;) =ek^{x). 
The mistuned version of the PDE ( fTSl l is then given by 



d'^v dv 2 



km dv kg d^v 
po dx 2pQ dx"^ 



(28) 



We study the problem of improving the stability margin by 
judicious choice of k„i{x) and ks{x). The results of our in- 
vestigation, carried out in the following sections, provide a 
systematic framework for designing control gains in the platoon 
by introducing small changes to the symmetric design. 
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A. Mistuning-based design for scenario I 

The control objective is to design mistuning profiles km{x) 
and ks{x) to minimize the least stable eigenvalue . To achieve 
this, we first obtain an explicit asymptotic formula for the 
eigenvalues when a small amount of asymmetry is introduced 
in the control gains (i.e., when e is small). For scenario I, the 
result is presented in the following theorem. The proof appears 
in Appendix II-BI 

Theorem 1: Consider the eigenvalue problem for the mis- 
tuned PDE ( |28] | with Dirichlet boundary condition STl\ corre- 
sponding to scenario I. The eigenvalue pair is given by the 
asymptotic formula 



Fig. 5. Mistuned front and back gains k^^^ and A;'''' of tlie vehicles in 
a platoon with fcp = 1 and e = 0.1. Figure (a) shows the gains chosen 
according to Corollary |2] to be optimal for scenario I for small e: k^^^ = 



ko (1 + 0.1(2_H"(7r - iS) - 1)) , k. 



ko (1 - 0.1(2_H"{7r - i<5) - 1)), 



where H{-) is the Heaviside function and S is defined in ^5)- Figure (b) shows 
the optimal mistuned gains for scenario II with the same parameters, which turns 



out to be (see Corollary |3) = l.lfco and fcj"' = 0.9feo for i = 1 



,Af. 



10" 



,-1 e 



'2boN J„ 



2ir 



km{x) sv[Y{lx)dx + O(e^) 



« * © 



I 



2n 



k„i{x) sm{lx)dx + O(e^) + 0(-^),C^ lo" 



'2boN Jo 

that is valid for each / in the limit as e ^ and N ^ oo. □ 
It is apparent from the Theorem above that to minimize 
the least stable eigenvalue s^, one needs to choose only k^ 
carefully; kg has only 0{-^) effect. Therefore we choose 
ks{x) = 0, or, equivalent^, fc(/^P"'-*)(x) = -^(''■p"''*) (a;), which 
leads to fcm(x) = 2k'^J''P^^*'\x). The most beneficial control 
gains are now can be readily obtained from Theorem [T] which 
is summarized in the next corollary. 

Corollary 2 (Mistuning profile for Scenario I): Consider the 
problem of minimizing the least-stable eigenvalue of the 
PDE (|28] ) with Dirichlet boundai-y condition (fTTl ) by 
choosing fe'^'^"''*) (x) e i°°([0, 27r]) with norm-constraint 
= max^e[02^] /c(/'P"''*)(x) = 1 and 
^(b,p«rt)(2;) = _k'^f,purt)f^^y jjjg jjjj^jj as e 0, the optimal 

mistuning profile is given by k^^ '"p^^^^x) — 2{H{x — -k) ~ i), 
where H{x) is the Heaviside function: H{x) — 1 for a; > 
and H{x) = for x < 0. With this profile, the least stable 
eigenvalue is given by the asymptotic formula 



hoN 



in the limit as e — > and iV — > oo. 



□ 



The result shows that even with an arbitrarily small amount of 
mistuning e, one can improve the closed-loop platoon damping 
by a large amount, especially for large values of N . The least- 
stable eigenvalue asymptotes to as 0{j^) in the mistuned 
case as opposed to 0{-^) in the symmetric case. 

Figure |5ja) shows the gains for the individual vehicles (that 
are obtained from sampling the functions fc''^^(x) and k'^^^x)), 
suggested by Corollary |2] for a 20 vehicle platoon, with ko = 1 
and e = 0.1: 



„(/) 



l + 0.2(i?(7r-i^) -0.5), and 



kl"' = l-0.2{H{n-iS)-0.5), 

where S is the desired inter-vehicular spacing in the scaled 
y coordinates, and is defined in (|5]l. A confirmation of the 
predictions of Corollary |2] is presented in Figure |6] Numerically 
obtained mistuned and nominal eigenvalues for both the PDE 
and the platoon state-space model are shown in the figure, with 
mistuned gains chosen as shown in Figure|5la). The figure shows 
that 



symmetric bidi. (platooif) 
symmetric bidi. (PDE) 
mistuned (platoon) 
mistuned (PDE) 
- Corollary © 



50 100 
N 



Fig. 6. Stabihty margin improvement by mistuning in Scenario I. The figure 
shows the least stable eigenvalue of the closed loop platoon (i.e., of A]^_p in 
(To)) and of the PDE (28) with Dirichlet boundary conditions, with and without 
mistuning, for a range of values of N. Parameters for the nominal case are 
ko = 1 and bo = 0.5, and the mistuning amplitude is e = 0.1. The mistuned 
control gains are shown in Figure [Sja)- The legend "Corollary |2]' refers to the 
prediction by Corollary |2] for large A^. 



1) the platoon eigenvalues match the PDE eigenvalues accu- 
rately over a range of N, and 

2) the mistuned eigenvalues show large improvement over 
the nominal case even though the controller gains differ 
from their nominal values only by ±10%. The improve- 
ment is particularly noticeable for large values of N, while 
being significant even for small values of N. 

For comparison, the figure also depicts the asymptotic eigen- 
value formula given in Corollary |2l 

Figure |4|b) graphically illustrates the mechanism by which 
mistuning affects the movement of eigenvalues sf eis N in- 
creases. By properly choosing the mistuning patterns k,n{x) and 
ks{x), damping can be "exchanged" between the eigenvalues 
and sj~ so that the less stable eigenvalue "gains" stability at 
the expense of the more stable eigenvalue . The net amount 
of damping is preserved, since + = —bo (as seen from 
Theorem [TJ. 



B. Mistuning-based design for scenario II 

For scenario II, asymptotic formula for the eigenvalue (coun- 
terpart of Theorem [Til is summarized in the following theorem. 
The proof is entirely analogous to the proof of Theorem [T] and 
is therefore omitted. 

Theorem 2: Consider the eigenvalue problem for the mis- 
tuned PDE ( |28] l with Neumann-Dirichlet boundary condi- 
tion (fTSl ) corresponding to scenario II. The eigenvalue pair 
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is given by the asymptotic formula 



I 



4fooA^ JO 
I 



(:r)sin(^)dx + 0(62) + 0(^), 



10" 



27r 



that is vaHd for each / in the limit as e 



(x) sm( — )dx 



and N ^ oo. 



□ 



As with scenario I, here again we use the above result to 
determine the most beneficial profile km{x) for small e: 

Corollary 3 (Mistuning profile for Scenario II): Consider 
the problem of minimizing the least-stable eigenvalue of the 
PDE (|28] ) with Neumann-Dirichlet boundary conditions ( fTsT i 
by choosing £ L°°([0,27r]) with norm-constramt 

max^g[o,2^] fcW'P"''*)(a;) = 1, and fc(^'P'"-*)(a;) = -fc(^'P'"'*)(a;). 
In the limit as e 0, the optimal k^S'P^^*') is given by 
k'^S 'P^^^) {x) — 1. With this profile, the least-stable eigenvalue 
is given by the asymptotic formula 



CO 
SO 

I 



in the limit as e 



and iV -> oo. □ 



The result shows that, as in scenario I, it is possible to 
improve the closed-loop stability margin in scenario II with an 
arbitrary small amount of mistuning e such that the least-stable 
eigenvalue s'l asymptotes to as 0[j^) in the mistuned case as 
opposed to 0{^) in the symmetric case. The gains suggested 
by Corollary [3] with fco = 1 and e = 0.1 are: 



fcK^ = 1.1, 



and 



0.9, 



which are shown in Figure |5jb). Numerically obtained least 
stable eigenvalues for the PDE and the platoon state-space model 
for scenario II are shown in Fig. |2]for a range of values of N . 
It is clear from the figure that, as in scenario I, the mistuned 
eigenvalues show an order of magnitude improvement over their 
values in the symmetric bidirectional case with only ±10% 
variation. 

Remark 2 (Robustness to small changes from the optimal gain 
An advantage of the mistuning design is that mistuned closed 
loop eigenvalues are robust to small local discrepancies in 
the control gains from the optimal ones. This can be seen 
(for scenario I) from the asymptotic eigenvalue formulas of 
Theorem[Tl which shows that one would obtain a 0(-^) estimate 
for any choice of k„i{x) such that J^^ km{x) sm{x)dx 7^ 0. A 
similar argument holds for scenario II. 



C. Simulations 

We now present results of a few simulations that show the 
time-domain improvements - manifested in faster decay of 
initial errors - with the mistuning-based design of control gains. 
Simulations were carried out for a platoon of = 20 vehicles 
with scenario I, i.e., with fictitious lead and follow vehicles. 
The desired gap was A = 1 and desired velocity was Vd ~ 5. 
The initial velocity of every vehicle was chosen as the desired 
velocity and the initial position of the vehicle was chosen as 
Zi{0) = iA - 0.5 for i = {1, . . . , N}. As a result, the initial 
relative position error and velocity error of every vehicle was 



symmetric bidi.°(p}atoon) 
symmetric bidi. (PDIi) 
mistuned (platoon) 
mistuned (PDE) 
-Corollary ^ 

20 50 1 

N 



Fig. 7. Stability margin improvement by mistuning in scenario II. The figure 
shows the least stable eigenvalue of the closed loop platoon (i.e., of Aj^ in \10\ ) 
and of the PDE )28t with Neumann-Dirichlet b.c, with and without mistuning, 
for a range of values of A^. The parameters for the nominal case are fco = 1 
and bo = 0-5, and the mistuning amplitude is e = 0.1. The mistuned control 
gains that are used are shown in Figure |5lb). The legend "Corollary [3]' refers 
to the prediction by Corollary |3] of mistuned PDE eigenvalues. 



zero except for the first vehicle, whose relative position error 
with respect to the fictitious lead vehicle was 0.5. 

Figure |8] shows the time-histories of the absolute and relative 
position errors of the individual vehicles with a symmetric 
bidirectional control, where the control gains were chosen as 
k^^ = fcf ^ = 1 and b, = 0.5 for i = {1, . . . , 20}. The absolute 
position error of the vehicle is Zi — Zf and the relative 
position error is Zi^i — Zi — A. 

Figure |9] shows the time-histories of the absolute and relative 
position errors for the platoon with mistuned controller gains. 
The mistuning gains used for the simulation are the ones 
shown in Figure [Sja) (chosen according to Corollary |2]i so that 
maximum and minimum gains over all vehicles is within ±10% 
of the nominal value. On comparing Figures [8] and |9l we see 
'^ihat the errors in the initial conditions are reduced faster in the 
mistuned case compared to the nominal case. These observations 
are consistent with the improvement in the closed-loop stability 
margin with the mistuned design. 

VI. Discussion on mistuning design 

There are several remarks to be made regarding the mistuning 
based design. We first comment on the implementation issues, 
in particular, on the effect of small platoon size on the proposed 
design, and on the information requirements for its implemen- 
tation. 



A. Large vs. small N 

The PDE model is developed for large N. However, detailed 
numerical comparison between the PDE and the discrete state 
space model shows that the PDE model provides quantitatively 
correct predictions even for small values of N (see Figures [5] |6] 
and |7]i. The PDE has an infinite number of eigenvalues as 
opposed to a finite number for the discrete platoon. So, one can 
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not expect an exact match. However, PDE eigenvalues exactly 
match the least stable and other dominant eigenvalues of the 
discrete platoon (see Figured and Figure [Tol). In a similar vein, 
the benefits of mistuning are also realized for small values of N. 
For example, when the number of vehicles is 20, a mistuning 
of ±10% results in an improvement in the stability margin - 
as measured by the real part of the least stable eigenvalue - 
of 150% (from -0.0491 to -0.1281 ) in scenario I and an 
improvement of 400% (from —0.012 to —0.05) in scenario II 
over the symmetric case. 

B. Information requirements 

In order to implement the beneficial mistuned controller gains 
designed above, every vehicle needs the following information 
(in addition to what is needed to use a symmetric bidirectional 
control): (1) the mistuning amplitude e, and (2) in scenario I, 
whether it is in the front half of the platoon or not. This 
information can be provided to the vehicles in advance. In 
scenario II, only the value of e is needed. 

It is possible that due to vehicles leaving and joining the 
platoon, information on whether a vehicle belongs to the front 
half of the platoon may become erroneous with time, especially 
for the vehicles that are close to the middle. In scenario I, such 
error may lead to a non-optimal gains used by the vehicles. 



However, since the improvement in closed loop stability margin 
due to mistuning is robust to small deviations in the gains from 
the optimal ones (see Remark |2]i, errors in determining whether 
a vehicle belongs to the front half of the platoon or not will not 
greatly affect the improvement in stability margin. Note that in 
scenario II this issue does not even arise. 



C. Large asymmetry 

Although the mistuning profiles described in Corollaries |2] 
and |3] are optimal in the limit as e 0, one would like to be 
able to use them with somewhat larger values of e to realize 
the benefit of mistuning. To do so, one has to preclude the 
possibility of "eigenvalue cross-over", i.e., of the second (s^) 
or some other marginally stable eigenvalue from becoming the 
least stable eigenvalue in the presence of mistuning. It turns 
out that such a cross-over is ruled out as a consequence of 
the Strum-Liouville (S-L) theory for the elliptic boundary value 
problems. The standard argument relies on the positivity of the 
eigenfunction corresponding to s^; the reader is referred to ||33l 
for the details. Figure [TO] verifies this numerically by depicting 
the six eigenvalues closest to (for both the PDE and the 
discrete platoon) as a function of N when mistuning is applied. 
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Fig. 10. The real parts of six eigenvalues (closest to 0) of the closed loop 
platoon dynamics for Scenario I, and their comparison with the PDE eigenvalues 
with Dirichlet-Dirichlet boundary conditions, with controller gains mistuned 
as those shown in Figure \5\ As predicted by the S-L theory, the least stable 
eigenvalue stays the least stable, although eigenvalues that are more stable merge 
with it as N increases. 



Fig. 1 1 . Hoo norm of the transfer function Gwe from disturbance w to spacing 
en'or e in j29K with and without mistuning, for scenario I. The mistuned gains 
used are shown in Figure |5la). Nor ms are computed using the Control Systems 
Toolbox in MATLAB©. 



D. Sensitivity to disturbance 

Automated platoons suffer from high sensitivity to external 
disturbances; which is referred to as "string instability" or 
"slinky-type effects" ||20]| . |[l], ifTSll . Here we provide numerical 
evidence that mistuning also helps in reducing the sensitivity to 
disturbances. 

When external disturbances are present, we model the dy- 
namics of vehicle ihy Zt = Ui + Wi, where Wi is the external 
disturbance acting on the vehicle. In the y coordinates, the 
vehicle dynamics become yi — ut + m;, where Wi := 2nWi/L. 
In scenario I, the state space model of the entire platoon 
becomes. 



V' = A 



L-F 



(29) 



e := \e^^^ 



where = [y-^jv"^]^, Al-f, w = [wi,'W2, ■ ■ ■ ,wnY , and 

JfhT 



^ , . . . , ] is a vector of front spacing errors e,- 

Vi-i - Vi- 

The Hoo norm of the transfer function G^je from the distur- 
bance w to the inter-vehicle spacing errors e is a measure of 
the closed loop's sensitivity to external disturbances ||7l, ifTSl . 
Figure [TT] shows a plot of the iJoo norm of G^e as a function 
of N , with and without mistuning. The mistuning profile used 
is the same as the one used for the eigenvalue trends reported 
in Figure |6] It is clear from the figure that ±10% mistuning 
results in large reduction of the iJoo norm of G^e- Although this 
reduction is more pronounced for large N, it is still significant 
for small N . In particular, for iV = 20, a 10% mistuning yields 
approximately 50% reduction in the Hoo norm (from 6.69 to 
3.38). 

Apart from the Hoo norm of Gde, there are other ways to 
measure sensitivity to disturbances. In [21 1, the transfer function 
from disturbance acting on the lead vehicle to spacing error on 
the vehicle is analyzed. Detailed analysis of the effect of 
mistuning on sensitivity to disturbances will be a subject of 
future work. 



VII. Conclusion 

We developed a PDE model that describes the closed loop 
dynamics of an A^-vehicle platoon with a decentralized bidirec- 
tional control architecture. Analysis of the PDE model revealed 
several important features of the problem. First, we showed that 
when every vehicle uses the same controller with constant gain 
that is independent of N (the so-called symmetric bidirectional 
architecture), the least stable eigenvalue of the closed loop 
decays to as 0{-^). Second, and more significantly, analysis 
of the PDE suggested a way to ameliorate the progressive loss 
of stability with increasing N, by introducing small amounts 
of "mistuning", i.e., by changing the controller gains from their 
nominal symmetric values. We proved that with arbitrary small 
amounts of mistuning, the decay of the least stable closed loop 
eigenvalue can be improved to 0(-^). Several comparisons with 
the numerically computed eigenvalues of state-space model of 
the platoon confirm the predictions of the PDE-based analysis. 

Although the PDE model is derived under the assumption that 
the number of vehicles, N, is large, in practice the PDE pro- 
vides quantitatively correct predictions for the discrete platoon 
dynamics even for relatively small values of N . The amount of 
information that is needed to implement the mistuned control 
gains (over that in the symmetric bidirectional architecture) is 
quite small and need to be provided only once. Furthermore, 
the stability improvement due to mistuning is robust to small 
errors (between the actual gains used and the optimal mistuned 
gains) that may occur in practice due to changes in the number 
of vehicles in the platoon over time. 

The advantage of the PDE formulation is reflected in the ease 
with which the closed loop eigenvalues are obtained for two 
different boundary conditions, with lead and follow vehicles as 
well as with only a lead vehicle. Certain important aspects of 
the problem, such as the beneficial nature of forward-backward 
asymmetry in control gains, is revealed by the PDE while 
they are difficult to see with the (spatially) discrete, state-space 
model. 

Numerical calculations show that the mistuning design also 
reduces sensitivity to disturbances of the closed-loop platoon. 
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Analysis of the beneficial effect of mistuning in reducing 
sensitivity to external disturbances is a subject of future research. 
In the future, we also plan to examine PDE-based models for 
modeling and analysis of fleet of vehicles as in 2 or 3 spatial 
dimensions. 

References 

[1] S. Darbha, J. K. Hedrick, C. C. Chien, and P. loannou, "A comparison of 
spacing and lieadway conti'ol laws for automatically controlled vehicles," 
Vehicle System Dynamics, vol. 23, pp. 597-625, 1994. 

[2] J. K. Hedrick, M. Tomizuka, and P. Varaiya, "Control issues in automated 
highway systems," IEEE Control Systems Magazine, vol. 14, pp. 21 - 32, 
December 1994. 

[3] R. E. Chandler, R. Herman, and E. W. Montroll, "Traffic dynamics: Studies 
in car following," Operations Research, vol. 6, no. 2, pp. 165-184, Mar. 
- Apr. 1958. 

[4] W. S. Levine and M. Athans, "On the optimal eiTor regulation of a string of 

moving vehicles," IEEE Transactions on Automatic Control, vol. AC- 11, 

no. 3, pp. 355-361, July 1966. 
[5] S. M. Melzer and B. C. Kuo, "A closed-form solution for the optimal 

error regulation of a string of moving vehicles," IEEE Transactions on 

Automatic Control, vol. AC- 16, no. 1, pp. 50-52, February 1971. 
[6] M. R. Jovanovic and B. Bamieh, "On the ill-posedness of certain vehicular 

platoon control problems," IEEE Transactions on Automatic Control, 

vol. 50, no. 9, pp. 1307 - 1321, September 2005. 
[7] P. Seller, A. Pant, and J. K. Hedrick, "Disturbance propagation in vehicle 

strings," IEEE Transactions on Automatic Control, vol. 49, pp. 1835-1841, 

October 2004. 

[8] J. D. Wolfe, D. F. Chichkat, and J. L. Speyer, "Decentralized controllers 
for unmaned aerial vehicle formation flight," in Guidance, Navigation and 
Control Conference, 1996, pp. July 29-31. 
[9] P. K. C. Wang, F. Y. Hadaegh, and K. Lau, "Synchronized formation 
rotation and attitude control of multiple free-flying spacecraft," Journal of 
Guidance, Control, and Dynamics, vol. 22, no. 1, pp. 28-35, 1999. 

[10] S. S. Stankovic, M. J. Stanojevic, and D. D. Siljak, "Decentralized 
overlapping control of a platoon of vehicles," IEEE Transactions on 
Control Systems Technology, vol. 8, pp. 816-832, September 2000. 

[11] P. Li and A. Shrivastava, "Traffic flow stability induced by constant 
time headway policy for adaptive cruise control vehicles," Transportation 
Research Part C: Emergent Technologies, vol. 10, pp. 275-301, 2002. 

[12] L. E. Peppard, "String stability of relative-motion FID vehicle control 
systems," IEEE Transactions on Automatic Control, pp. 579-581, October 
1974. 

[13] P. Barooah and J. P. Hespanha, "Error amplification and distrubance 
propagation in vehicle strings," in Proceedings of the 44th IEEE conference 
on Decision and Control, December 2005. 

[14] K. C. Chu, "Decentralized control of high-speed vehicle strings," Trans- 
portation Science, vol. 8, pp. 361-383, 1974. 

[15] Y. Zhang, E. B. Kosmatopoulos, P. A. loannou, and C. C. Chien, 
"Autonomous intelligent cruise control using front and back information 
for tight vehicle following maneuvers," IEEE Transactions on Vehicular 
Technology, vol. 48, pp. 319-328, January 1999. 

[16] S. E. Shladover, "Longitudinal control of automotive vehicles in close- 
formation platoons," Journal of Dynamic systems, Measurements and 
Control, vol. 113, pp. 302-310, December 1978. 

[17] H.-S. Tan, R. Rajamani, and W.-B. Zhang, "Demonstration of an automated 
highway platoon system," in American Control Conference, vol. 3, June 
1998, pp. 1823 - 1827. 

[18] X. Liu, S. S. Mahal, A. Goldsmith, and J. K. Hedrick, "Effects of 
communication delay on string stability in vehicle platoons," in IEEE 
International Conference on Intelligent Transportation Systems (ITSC), 
August 2001. 

[19] P. Barooah, P. G. Mehta, and J. P. Hespanha, "Control of large vehicular 

platoons: Improving closed loop stability by mistuning," in The 2007 

American Control Conference, July, pp. 4666-4671. 
[20] S. Darbha and J. K. Hedrick, "String stability of interconnected systems," 

IEEE Transactions on Automatic Control, vol. 41, no. 3, pp. 349-356, 

March 1996. 

[21] R. H. Middleton and J. H. Braslavsky, "String instability in 
classes of linear time invariant formation control with Hmited 
communication range," 2008, submitted for pubHcation. [Online]. 
Available: http://www.hamilton.ie/rick/publications/StringStability.pdf 

[22] S. K. Yadlapalli, S. Darbha, and K. R. Rajagopal, "Information flow and 
its relation to stability of the motion of vehicles in a rigid formation," 
IEEE Transactions on Automatic Control, vol. 51, no. 8, August 2006. 



[23] B. Shapiro, "A symmetry approach to extension of flutter boundaries via 
mistuning," Journal of Propulsion and Power, vol. 14, no. 3, pp. 354-366, 
1998. 

[24] O. O. Bendiksen, "Localization phenomena in structural dynamics," 

Chaos, Solitons, and Fractals, vol. 11, pp. 1621-1660, 2000. 
[25] A. J. Rivas-Guerra and M. P. Mignolet, "Local/global effects of mistuning 

on the forced response of bladed disks," Journal of Engineering for Gas 

Turbines and Power, vol. 125, pp. l-Il, 2003. 
[26] P. G. Mehta, G. Hagen, and A. Banaszuk, "Symmetry and symmetry 

breaking for a wave equation with feedback," SIAM Journal of Dynamical 

Systems, vol. 6, no. 3, pp. 549-575, 2007. 
[27] M. Lighthill and G. Whitham, "On kinematic waves II: a theory of traffic 

flow on long crowded roads," in Royal Society, London Series A, 1955. 
[28] D. Helbing, "Traffic and related self-driven many-particle systems," Review 

of Modern Physics, vol. 73, pp. 1067-1141, 2001. 
[29] D. Jacquet, C. C. de Wit, and D. Koenig, "Traffic control and monitoring 

with a macroscopic model in the presence of strong congestion waves," 

in 44th IEEE Conference on Decision and Control & European Control 

Conference, 2005, pp. 2164-2169. 
[30] P. Y. Li, R. Horowitz, L. Alvarez, J. Frankel, and A. M. Robertson, 

"An automated highway system link layer controller for traffic flow 

stabilization," Transportation Research, Part C, vol. 5, no. 1, pp. 11-37, 

1997. 

[31] L. Alvarez, R. Horowitz, and P. Li, "Traffic flow control in automated 
highway systems," Control Engineering Practice, vol. 7, pp. 1071-1078, 
1999. 

[32] C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral 

Methods in Fluid Dynamics, ser Springer Series in Computational Physics. 

New York: Springer- Verlag, 1983. 
[33] L. C. Evans, Partial Differential Equations, ser. Graduate Studies in 

Mathematics. American Mathematical Society, 1998, vol. 19. 
[34] A. Pazy, Semigroups of linear operators and applications to partial 

differential equations, ser Applied Mathematical Sciences. New York: 

Springer-Verlag, 1983, vol. 44. 

Appendix I 
Technical results 

A. Solution properties of PDE il5\ . 

In this section, we use the semigroup theory to obtain results 
on well-posedness of the PDE ( fTSl l. To apply these methods, we 
first re-write the PDE as a first order evolution equation: 

dp _ _ dv_ 

dt ~ , POox 



dt 



bv 



A 



P 

V 

"(30) 

where A is a linear operator; ko{x) :— fc+(x) and ki{x) :— 
k~{x) — 2^^i^(a^)- We will assume these coefficients 
ko{x),ki{x) £ i°°([0, 27r]) and ko{x) > 0. p has the units 
of and the physical interpretation of density perturbation. 
Using dSOl l. we denote the initial/boundary value problem as: 

z{x,t) = Az{x,t) foTxeX, t>0 

z{x,0) = ^0(2;), (31) 

where z{x,t) :— [p{x,t),v{x,t)], zo{x) = [po{x),vo{x)] and 
A is defined in dSOl i: po and vq will be assumed to functions 
in appropriately defined Banach spaces. The main goal of this 
section will be to show that the solution for the linear prob- 
lem ( l30b can be expressed in terms of a semigroup provided 
eigenvalues of the operator A satisfy appropriate bounds. We 
begin with a discussion of the notation. 

Preliminaries and Notation. We denote z :== [p,v], L'^{X) 
denotes the Hilbert space of square integrable functions on X 
(II^IIl^ / "^"^dx), denotes the Sobolev space of functions 
such that derivatives up to fc*'*-order exist in a weak sense and 
belong to L'^{X) (the Sobolev norm is denoted by || • ||jjfc), and 
Hq denotes the Sobolev space of functions that satisfy the 
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Dirichlet boundary condition. We denote Z := x L^, and 
equip it with a norm || • ||. Let V{A) x {H^ n L^) and 

consider the right hand side of evolution equation dSOl l as an 
unbounded but closed densely defined linear operator 



dv 



A : V{A) C Z ^ Z. 



(32) 



A real number s belongs to p{A), the resolvent set for A, 
provided the operator si — A : I? (A) Z is 1-1 and onto. 
For s e p{A), the resolvent operator Rg := {si — A)^^. 
Finally, we recall that a one-parameter family of linear operators 
{S{t)}t>a is a C°-semigroup if 1) S{0)z = z for all 2 e 2) 
S{t + s)z = S{t)S{s)z for all i, s > and z e Z, and 3) the 
mapping t S{t)z is continuous from [0, 00) into Z. A C'^ 
semigroup is a contraction semigroup if || 5(^)2:11 < ||z|| for all 
t > 0. The Hille-Yosida theorem states that a closed densely 
defined linear operator A is the generator of a contraction 
semigroup if and only if 



(0,00) Cp(A) and p^zll < Vz G Z. 

s 



(33) 



Our strategy will be to apply Hille-Yosida theorem to deduce 
solution properties of the evolution equation ( l3Tl i. Following 
closely the development in |.33„|. there are three steps to accom- 
plish this: 1) we show that ^ is a densely defined closed linear 
operator on Z, 2) characterize the resolvent set by considering 
the eigenvalue problem, and 3) show the bound ( [33] l for the 
resolvent. Step 2 will lead to an eigenvalue problem, whose 
analysis and optimization is the subject of this paper We present 
details for the three steps next: 

1) The domain of A, T){A), is dense in Z because is 
dense in L^. To show A is closed, consider a sequence 
{P7n,v„i} C V{A) such that 

iPm,Vm) ip,v) (34) 

A{p^,v„r) {f,g), (35) 

where the arrow notation denotes the fact that the conver- 



gence is in Z ^ L X L . Since v„ 



V so -po^ = 



f £ I?, i.e., V e H^. Now, is Cauchy in I? by 
and {^} is Cauchy in L"^ by (O and 

dvm _ dm 
dx 



dx 



UU[ .. .. 



2) 



(36) 

so {vm} is Cauchy in and w,„ ~> v. By repeating 
essentially the same argument, one also finds that p E 

and pm P- Consequently, A{pm,Vm) ^ A{p,v) and 
A{p,v) = {f,g). 

Let s > 0, {f,g) e Z — L'^ X L^, and consider the 
operator equation 



p 




■ / " 


V 




. g _ 



(si - A) 

This is equivalent to two scalar equations 

dv 



Using the first equation to write sp = —po-g^ 
implies 

s'^v + bsv + Lv = h, 



where 



3) 



= / 



(37) 
(38) 



\kT,{x)p + -^^(fco(x)p) 
Po 2pg dx 



bv 



1 d 



dv , 

^di' 



1 

Pa 



dv , 
^di' 



/, this 
(40) 

(41) 



is an elliptic operator (because kQ{x) > for all x E X) 
and h ^ sg - -^^{ko{x)f) - jiki{x)f (note that 
h e II^^{X)). Consequently, solutions of (dJTb) can be 
studied in terms of solutions of ((|4TI)). The spectrum of A 
is completely characterized by the spectrum of L. We will 
obtain spectral bounds, dependent upon A:o(x) and ki{x), 
in the following sections. In particular, we will establish 
that Real[s] < a for some a < and thus p{A) D {a, 00). 
For ki{x) = 0, its turns out that [0, 00) C p{A) for 
any choice of positive kQ{x) (this is also clear from the 
symmetric eigenvalue problem (|4TI) ). 
If a positive s E p{A), there exists a unique solution 
(p, v) E Z for (|38])-(|391) via the theory of elliptic opera- 
tors: solve (|40ll to obtain v E Hq and sp = -poff +/■ We 
write the solution as {p,v) ~ Rs{f,g), define a bilinear 
form 

1 



2p; 



B[p,s] := — ^ / ko{x)p{x)s{x)dx, 



(42) 



'0 Jx 



for p,s E and consider an equivalent norm (on Z) for 
solutions {p, v) as: 



\\{p,v)\\ := B[p,p] + \\v\\l2 



(43) 



To obtain the resolvent bound, we multiply ( [39l ) by v and 
use integration by parts: 



B[p, p]) + b\\v\\L2 + — J ki{x)pvdx 

= j gvdx + B[pJ]. 

In general, the bound depends upon ki{x). For ki{x) = 0, 
we have 

smv)\\' <{s + b)\\v\\L2+B[p,p]) = 

gvdx + B[pJ]<\\{f,g)\\-\\{p,v)\\, 

where the first inequality holds because s > and 
6 > and the last inequality follows from the generalized 
Cauchy-Schwarz inequality. As a result, ||i?s(/, < 

and ||i?.||<i. 
For the general case where ki{x) is not identically zero, 
one expresses the operator 



= g (vEL^nH'o)- 
(39) 



where 

Ao 

A 







p 




V 




p 




v 





A = Ao + A, 



(44) 



-POd^ 

-b 









p 




V 



-jikiix) 





p 




v 
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In words, Aq is the operator with ki{x) ~ and A is 
the operator due to ki{x). We note that ^ is a bounded 
perturbation of Aq (on Z). We have already showed the 
existence of a C''-semigroup for Aq. For the general op- 
erator A, the existence follows from using a perturbation 
theorem (see Theorem 1.1 in Ch. 3 of Ii34l ). 

B. Proof of Theorem Q] 

Proof of Theorem [7] The spatial inhomogeneity introduced by 
the x-dependent coefficients k„i{x) and ks{x) destroy the spatial 
invariance of the nominal PDE ( |20] |. Hence, the Fourier basis 
- eigenfunctions of the Laplacian - no longer lead to a diag- 
onalization of the mistuned PDE. The methods of section |IV] 
thus need to be suitably modified. In order to compute the 
eigenvalues for the mistuned PDE ( |28] ), we take a Laplace 
transform of ( l28T l and get 



,5^ 



s^i] + bosrj = e 



km dr] 

Po dx 



ks d'^-q 



(45) 



where ri{x) is the Laplace transform (with respect to t) of 
v{x,t). We are interested in eigenvalues of (05]) with Dirichlet 
boundary conditions, i.e., the values of s for which a solution 
to the homogeneous PDE i45[ exists with boundary conditions 
77(0) = ?7(27r) = 0. To obtain these eigenvalues, we use a per- 
turbation method expressing the eigenfunction and eigenvalue 
in a series form: 



?7(x) = ?7o(a;) + e?7i(a;) + O(e^), 



ro + eri+0{e^). 



(46) 



We note that eri denotes the perturbation to the nominal 
eigenvalue tq as a result of the mistuning. Substituting ( |46] | 
in (|45] | and doing an 0(1) balance, we get 



0(1): -al{rio)xx+rlr]o + brorio = 0, 
whose eigen-solution is given by 



Vo 



di sin 



where < •, • > denotes the standard inner product in ^^(o, 2tt). 
This leads to an equation 



(2ro + bo)ri 



I 



47rpo Jo 



km{x) sm{lx)dx 



f^'^ Ix 
2 / ks{x)siu'{-)dx (49) 
'0 Jo ^ 



For values of To — sf(fS), where sfifS) is given by ( |24] |. the 
equation above leads to an expression for perturbation in the 
two eigenvalues. We denote these perturbations as . For tq = 
s^(0), we have from from Lemma [T] that ho >> |2ro| when 
I « Ic, which happens for every Z as — > cx) (see eq. dZSt). 
so that 



47rp( 



— / fc„(.T)sin(Zx)dx + 0(^). (50) 

Note that we have dropped the second integral on the right hand 
side of (|49]l because = 0{l/N'^) for large N. For tq = 

s^(0), 2ro 



Pa 

-2ho for I «lc and 

r2iT 



Note that 



— / km{x)iiin{lx)dx + 0{—^). 



= 0. 



(51) 



Putting the formulas for the perturbation to the eigenvalues 
and dsTl l in ( |46] |. we get 

s+ie) « s+(0) + ^ km{x) smilx)dx + 0{e') + O(^), 

I f^'^ 1 
s;-(e) w -60 - £ ^^^^^^ y fcm(a;) sin(/x)dx + ©(e^) + 0(-^). 

Since s'^{0) = 0{^) for Z < (Lemma [TJ and po = the 



(47) result follows. 



where I — 1,2, di is an arbitrary real constant, and S;^(0) 
is given by ( l24l i. Next, 



0(e): 



'9x2 



+ (^-o + ^o^-o) m 



km drjo 
Po dx 
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+ ^^^^0 ~ (2rori + 6ori)77o =: R 

Substituting ro = S/^(0) on the left hand side leads to a 
resonance condition for the right hand side term, denoted by 
R. In particular for a solution 771 to exist, R must lie in the 
range space of the linear operator 



92 



(48) 



For this self-adjoint operator, the range space is the complement 
of its null space {sin(Y)}. This gives the resonance condition 



as 



(^,sin(^)) 



